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Symmetry protected topological (SPT) phases with unusual edge excitations can emerge in 
strongly interacting bosonic systems and are classified in terms of the cohomology of their sym- 
metry groups. Here we provide a physical picture that leads to an intuitive understanding and 
wavefunctions for several SPT phases in d — 1, 2, 3 dimensions. We consider symmetries which 
include a Z2 subgroup, that allows us to define domain walls. While the usual disordered phase is 
obtained by proliferating domain walls, we show that SPT phases are realized when these prolif- 
erated domain walls are 'decorated', i.e. are themselves SPT phases in one lower dimension. For 
example a d — 2 SPT phase with Z2 and time reversal symmetry is realized when the domain walls 
that proliferate are themselves in a d — 1 Haldane/AKLT state. Similarly, d — 3 SPT phases with 
Z2 X Z2 symmetry emerges when domain walls in a = 2 SPT with Z2 symmetry are proliferated. 
The resulting ground states are shown to be equivalent to that obtained from group cohomology and 
field theoretical techniques. The result of gauging the Z2 symmetry in these phases is also discussed. 
An extension of this construction where time reversal plays the role of Z2 symmetry allows for a 
discussion of several d = 3 SPT phases. This construction also leads to a new perspective on some 
well known d = 1 SPT phases, from which exactly soluble parent Hamiltonians may be derived. 
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I. INTRODUCTION 

Symmetry protected topological (SPT) phases are 
gapped quantum phases with topological properties pro- 
tectected by symmetry^. The ground states of SPT 
phases contain only short-range entanglement and can 
be smoothly deformed into a totally trivial product state 
if the symmetry requirement is not enforced in the sys- 
tem. However, with symmetry, the nontrivial SPT order 
is manifested in the existence of gapless edge states on 
the boundary of the system which cannot be removed as 
long as symmetry is not broken. Many SPT phases have 
been discovered over the past decades^^^ and a general 
structure of the theory of SPT phases is emerging. In Id, 
SPT phases have been completely classified for general 
interacting bosonic/fermionic systems which carry non- 
trivial projective representations of the symmetry in their 
degenerate edge states^^^^. In fact, the Haldane/AKLT 
phase of spin-one quantum antiferromagnet, is a physical 
realization of a ID SPT phase, which is protected by spin 
rotation or time reversal symmetries. Two dimensional 
SPT phases were first discovered in topological insulators 
and superconductors which has subsequently been gen- 
eralized to three dimensions^^^^. Such free fermion SPT 
phases have realized experimentally^^^^ and also classi- 
fied completely^^'^"^. Recently, it was realized that two 
and higher dimensional SPT phases exist not only in free 
fermion systems but also in strongly interacting boson 
systems and a systematic construction is given based on 
the group cohomology of the symmetry^. While the con- 
struction provides a fixed point description in the bulk, 
it is hard to access edge dynamics and hence how the 
system responds to physical perturbations. Also it is 
not clear in what physically realistic systems can these 
bosonic SPT phases be realized. 



Much progress has been achieved recently in under- 
standing the low energy physics and finding physical re- 
alizations for some of the bosonic SPT phases. In 2D, 
a general understanding of SPT phases with Abelian 
(and time reversal) symmetry has been given in terms 
of Chern-Simons K-matrix theory^, which also provides 
a field theory of the protected edge states. A physi- 
cal realization of the U{1) SPT phases (which are ex- 
pected to have even integer quantized quantum Hall 
conductance^) has been proposed in bosonic cold atom 
systems with artificial gauge fields^. The 2D SPT phases 
with nonabelian S0{3) and SU{2) symmetry was stud- 
ied with nonlinear sigma model with quantized topo- 
logical terms and are found to have quantized spin 
transport^^. More recently, some 3D SPT phases have 
been understood within a field theoretic approach, which 
predicts surface vortices with projective representations 
and quantized magnetoelectric responses^. Recently, 3D 
SPT phases have been discussed from number of dif- 
ferent theoretical perspectives , including twisted vor- 
tex condensates^^ and the statistical magneto-electric 
effect^^. A useful perspective on SPT phases appears on 
gauging the symmetry discussed by Levin and Gu^ in 2D, 
and recently extended in 20, 21, 35-37 to 3D phases. Fi- 
nally, we note a special feature of 3D SPT phases is that 
their surface states could be gapped and fully symmetric 
if they develop topological order just at the surface^'^^^^. 

In this paper, we focus on SPT phases with symmetry 
group Z2 X G and present a simple construction of d di- 
mensions SPT phases by decorating the domain walls of 
Z2 configurations in the bulk with d—1 dimensional SPT 
states with G symmetry. Such a construction naturally 
reveals some of the special topological features of such 
phases. For example, it is easy to see that when a do- 
main wall is cut open at the boundary of the system, the 
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end points/loops of the domain wall carry gapless edge 
states of the d — 1 dimensional SPT state with G sym- 
metry. The same is true on the flux point/loops in the 
system when the Z2 part of the symmetry is gauged. In 
particular, we are going to present the construction of a 
2D SPT phase with Z2 x symmetry by decorating the 
ID Z2 domain walls with Haldane chains and a 3D SPT 
phase with Z2 x Z2 symmetry by decorating the 2D Z2 
domain walls with the nontrivial 2D SPT phase with Z2 
symmetry. These constructions hence provide a simple 
understanding of some higher dimensional SPT phases in 
terms of lower dimensional ones (which are usually better 
understood) and belongs to a hierarchical construction of 
SPT phases related to a hierarchy structure of group co- 
homology. These results are summarized in the figures 1, 
4, 8. Besides the Z2 x G symmetry, similar domain wall 
picture also allows us to construct SPT phases in 3D 
with Zj, Z2 X U{1) and x S0{3) symmetry, which 
produced a Z2 class in each case. 

The paper is organized as follows: In section II we 
present the construction for the 2D SPT phase with 
Z2 X Z2 symmetry while section III is devoted to the 
3D SPT model with Z2 x Z2 symmetry. We describe 
in detail for both ground state wavefunctions on a lat- 
tice, the symmetry action on the edge (and its relation 
to group cohomology) , and an low energy effective theory 
description of both the bulk and the edge which captures 
the special domain wall feature of the state. In section IV 
we discuss the consequence of gauging the Z2 symmetry. 
In section V and VI we construct ID SPT phase with 
Z2 X Z2 symmetry and some 3D SPT phases with time 
reversal symmetry respectively. The relation between the 
domain wall construction and the Kiinneth formula for 
group cohomology is discussed in section VII. We discuss 
the implications and the generalizations of our results in 
section VIII. Brief reviews of the group cohomoloy de- 
scription of SPT phases, the field theory description of 
2D and 3D SPT phases are given in the appendices. 



II. 2D SPT PHASE WITH Z2 x SYMMETRY 

Let's start with the simplest example in this construc- 
tion: a 2D SPT phase with Z2 x symmetry where Zj 
represents time reversal symmetry. Intuitively, the state 
is constructed by attaching ID Haldane chains with time 
reversal symmetry to the domain walls between Z2 con- 
figurations on the plane and then allow all kinds of fiuc- 
tuations in the domain wall configurations. With such a 
structure, it is then easy to see that when the system has 
a boundary where domain walls can end, the end point 
would carry a spin 1/2 degree of freedom (edge state of 
Haldane chain) which transform project ively under time 
reversal. 

In section II A, we describe a 2D lattice wave function 
of the ground state which is gapped and does not break 
any symmetry. Section II B then discusses what hap- 
pens on the boundary of the system. By identifying the 




FIG. 1. SPT phase ind — 2 with Z2 and time reversal symme- 
try (Z2 X Zj). A snapshot of the ground state wave function, 
where blue and grey are oppositely directed domains of the 
Z2 symmetry. The ground state preserves the Z2 symmetry 
since it is a superposition of domain configurations. The do- 
main walls themselves (black lines) are in a d = 1 SPT phase 
(the Haldane/ AKLT phase) protected by time reversal sym- 
metry. When they end at the edge of the system they create 
Kramers doublets, leading to a gapless edge state. 

relation between symmetry action on the edge and the 
nontrivial cocycle of Z2 x Zj symmetry group, we estab- 
lishes the existence of nontrivial SPT order in this model. 
It is known that this SPT phase with Z2 x Zj symmetry 
can be described with U{1) x U{1) Chern-Simons the- 
ory with a nontrivial symmetry action. We review this 
description in section II C and demonstrate the nontriv- 
ial topological feature of domain walls in the field theory 
language. 

A. Bulk Wave function on 2D lattice 




FIG. 2. Two-dimensional SPT model with Z2 x Zj symmetry, 
(a) each plaquette hosts a Z2 variable (big black dot) (b) each 
vertex hosts four spin 1/2 's (small blue dots) (c) two spin 1/2 's 
on the same link form a singlet if they are on a Z2 domain 
wall, other spin 1/2's form singlets within each vertex. 

Consider a honeycomb lattice as show in Fig. 2 (a) 
where each plaquette hosts a Z2 variable (big black dot) 
in state |0) or |1). The Z2 part of the symmetry fiips |0) 
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into |1) and |1) into |0). At each vertex, there are four 
spin 1/2's (one at the center and three on hnks as shown 
in Fig. 2 (b)). On each spin 1/2 time reversal symmetry 
acts as T = icFyK and satisfies = —1. On the total 
Hilbert space of each vertex, time reversal still satisfies 
= 1 and forms a linear representation of Zj. Now 
consider a Hamiltonian term V which enforces that in 
the ground state two spin 1/2's on the same link form 
a singlet pair if the link is on a Z2 domain wall while 
the remaining spin 1/2's which are not on a domain wall 
form singlets within each vertex (there are always an even 
number of these at each vertex). The effect of this term 
can be thought of as attaching Haldane chains to all the 
Z2 domain walls. A possible configuration in the ground 
state is shown in Fig. 2 (c) where the dotted lines rep- 
resent singlet pairing. A tunneling term between the Z2 
configurations ^ • is then added to the Hamiltonian 
which in the low energy sector of V flips Z2 variables to- 
gether with the related singlet configurations. Therefore, 
the ground state of the total Hamiltonian 



H = uV 



(1) 



when u is very large is an equal weight superposition 
of all possible consistent configurations of Z2 variables 
and singlets. The ground state is unique, gapped and 
preserves Z2 x symmetry. 



B. Edge state 



(a) 




and add a Hamiltonian term for them to be equal. This 
does not violate the symmetry of the system and is there- 
fore allowed.) The ID boundary degrees of freedom then 
contain two parts: the bond Z2 variables which come 
from plaquettes on the cut and the vertex spin variables 
when the neighboring bonds contain different Z2 config- 
urations, as shown in Fig. 3 (a). Symmetry acts by 
fiipping the Z2 variable and time reversing the spins. 

To study the boundary as an effective one dimensional 
system, we will use a slightly different description of the 
boundary degrees of freedom. In the previous descrip- 
tion, the existence of vertex degrees of freedom are de- 
pendent on the bond degrees of freedom. Equivalently, 
we can also think of the boundary as a local ID system 
with independent degrees of freedom on the bonds and 
the vertices, as shown in Fig. 3 (b). That is, we can imag- 
ine that there exist pseudo-spins (empty circles) on the 
vertices when neighboring Z2 variables are the same. On 
the pseudo-spins time reversal acts simply as T = cFxK 
and satisfies = 1. The total symmetry action is then 

Z2 : Tx on each r spin 



icjyK if (7 is on r domain wall 
(Ja^i^ot her wise 



It is easy to see that such a description of the bound- 
ary can be reduced to the previous one by polarizing 
each pseudo-spin in the z direction. Because = 1 on 
each pseudo-spin, polarizing them does not break time 
reversal symmetry. Using this description, we can see 
explicitly how this model is related to the nontrivial 3rd 
cocycle of the symmetry group and therefore contains 
nontrivial SPT order. We can also write down effective 
Hamiltonians for the boundary and solve for the low en- 
ergy dynamics, as we demonstrate in the following. 

In terms of the local degrees of freedom on the ID edge 
state, (J and r, the Z2 symmetry action on the edge reads 



n 



(2) 



FIG. 3. Two equivalent descriptions of the boundary of 2D 
SPT model with Z2 x symmetry, (a) Thick bonds rep- 
resent Z2 variable in state |1) and thin bonds |0), spin 1/2 
exists on their domain wall, (b) One Z2 variable per bond 
and one spin variable per vertex. Time reversal acts on spins 
in a way dependent of neighboring Z2 configurations. Dotted 
boxes represent local degrees of freedom on the ID boundary 
labeled by group elements. 

Interesting things happen when the system has a 
boundary. When the system is cut open, the Z2 do- 
main walls will have end points on the boundary. Be- 
cause the spin singlets are tied to domain walls in the 
bulk, on the boundary there will be isolated spin 1/2's 
on the end points of Z2 domain walls. Imagine that we 
cut the system through the middle of plaquettes. (In or- 
der not to cut through degrees of freedom, we can first 
split the Z2 variable into two, one on each side of the cut 



and time reversal acts as 



n(^ 



■TlT. 



i+1 



-CF, 



i+1 



r r 



^i<7;+i)i^ (3) 



We can write down an effective Hamiltonian satisfying 
the symmetry and solve for the dynamics on the ID edge. 
Some simple interaction terms which satisfy both the Z2 
symmetry (Eq.2) and the time reversal symmetry (Eq.3) 
include + cr*r^cr*+-^ and + r*~-^cr^r*. Therefore, a 
possible form of the dynamics of the edge is given by 
Hamiltonian 



(4) 



This Hamiltonian can be mapped exactly to an XY 
model by applying unitary transformations U = 



^-2 I 
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to each pair of 
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Ti and (7i variables. The Hamiltonian that we get after 
these transformations is 



H, = UHM- 



and the symmetry transformations are mapped to 



(5) 



(6) 



and 



Zl : n (^^^^4 + ^-^^^ial)K (7) 

i 

In order to see how the symmetry acts on the low en- 
ergy effective theory, we diagonalize the XY Hamilto- 
nian i^e, identify the free boson modes in the low en- 
ergy eigenstates and calculate the action of the sym- 
metry on the low energy states. The low energy 
states are labeled by quantum numbers Uk ^ Z and 
fik ^ k = 0,1,2..., where no and no label the 
total angular momentum and the winding number of 
the boson field respectively and Uk and fik^ k > 
0, label the left /right moving boson modes. The 
Z2 symmetry acts by mapping state |no, no, {n/^, n/c}) 
to (— )^fc>o |— no, —no, {n/e,n/e}) and time rever- 
sal symmetry acts by mapping state |no, no, {n/c, n/e}) 



to (-) 



no, — no, {n/e, n/c}), which is consistent 



with the low energy description given in the next section 
where Z2 symmetry acts as — (/>i, (/>2 ^ —4^2 and 

time reversal acts as + tt, (/)2 ^ (/>2 + tt on the 

chiral boson fields {^1,2}- 



C. Effective field theory description 

As discussed in Ref. 6, SPT phases in 2D may be 
characterized by the unusual transformation properties 
of their edge states, under the action of symmetry. The 
simplest situation, which captures a large fraction of SPT 
phases, is a Luttinger liquid edge with a single gapless 
bosonic mode. The edge fields are characterized by com- 
pact conjugate bosonic fields ^1, ^2 such that 



^i(x). 



27T 



— id{x — x') 



(8) 



Physically, e 



inserts a boson at point x along the 



edge, while e^^'^ inserts a phase slip. From the discus- 
sion in the previous section and in Appendix D, we see 
that the 2D SPT phase of interest implies the following 
transformation law: 



Z>1 ^ (/>1 + TT 
i)2 02 

/>2 ^ 02 + TT 



(9) 



Note, this edge cannot be gapped out without break- 
ing symmetry. Terms that gap the edge such as cos0i,2 
are forbidden by symmetry. If either one of the symme- 
tries is broken, then a trivial edge is possible. Finally let 
us verify that a domain wall of the Z2 order, terminat- 
ing at the edge carries a Kramer's pair, as the pictorial 
description implies. 

The operator that creates a domain wall, rotates 0i 
by TT everywhere (say) to the right of the point x. By 
the commutation relation 8, this is identified with the 
operator 



D{x) 



(10) 



Now, let us consider how the domain wall insertion op- 
erator transforms under time reversal T . In particular 
we would like to calculate the action of acting twice 

with time reversal and see if — ±1. Now D{x) 
g-*(02+7r)/2 ^ -iD''{x). Applying this again we find: 



D{x) A= -D{x) 



(11) 



thus we find that for this operator = —1. Hence it 
must carry a Kramer's degeneracy as expected from the 
physical picture. 



D. Connection to Group Cohomology 

Now let us establish the nontrivial SPT order in this 
state by identifying the connection of the symmetry ac- 
tions on the boundary with nontrivial third cocycles. 
When we think of the boundary as a local ID system 
with independent bond variables and vertex variables 
(i.e. treating spins and pseudo-spins as equivalent vari- 
ables), symmetry no longer acts on the degrees of freedom 
independently. The Z2 part of the symmetry Z2: Yli^x 
still acts on each Z2 variable independently. However, 
time reversal symmetry on the vertices depends on the 
Z2 configurations on the bonds, as shown in Eq.3. This 
is exactly the signature of SPT phases. 

Indeed, as was shown in Ref. 1, the boundary of dD 
STP phases can be thought of as a — ID local system 
where the symmetry acts in a non-onsite way. Without 
symmetry, the boundary can be easily gapped out. If the 
symmetry acts in an on-site way, the boundary can be 
simply gapped out by satisfying the symmetry on each 
site. However, with non-onsite symmetry related to non- 
trivial group cocycles, the boundary must remain gapless 
as long as symmetry is not broken. 

In particular, in the group cohomology construction as 
reviewed in appendix B, the boundary local degrees of 
freedom are labeled by group elements ai of the sym- 
metry group and the action of the symmetry operator 
involves two parts: 

• 1. changing the local states \ai) to \aai)^ 
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• 2. multiplying a phase factor given by nontrivial 
cocycles to each pair of ai and Q^i+i. 

More specifically, on the ID boundary of 2D SPT phases, 
the symmetry acts as 

0(a)|ai, ....,aAr) = ]^ /"(o^i, Q^i+i)|Q^Q^i, ...,aaAr) (12) 

i 

where is a phase factor given by the non- 

trivial 3rd cocycle ujs 

= us{a~^ai^i,a~l^a~'^,a) (13) 

Now we can show that the symmetry action on the 
boundary of the Z2 x model we constructed is of ex- 
actly this form. We can consider the pair of Z2 variables 
o^i = (dotted box in Fig. 3 (b)) as labeling group 

elements in the group Z2 x Zj. That is, we consider 
the Z2 state |0)/|1) of r as labeling the trivial/no ntrivial 
element in group Z2 and the spin state |t)/|i) of ^ la- 
beling the trivial/nontrivial element in group Zj. Then 
the Z2 part of the symmetry is the following mapping 

\Ti) ^ \fi) (14) 

and the time reversal part of the symmetry (given by 
either T = (JxK on the pseudo-spin or T = icFyK on the 
spin) also involves the mapping 

Wi) ^ \<yi) (15) 

Therefore, a general symmetry operation labeled by a = 
{t^(t)^ T ^ Z2^ (J ^ Zj, will first change group elements 
labels in each box 

Wi = {n^cFi)) \aai = (TTi.aai)). (16) 

Moreover, the time reversal symmetry also adds a (—1) 
phase factor when the neighboring Z2 variables are differ- 
ent and when the vertex variable was originally in state 
It)- Such a phase factor on spin can be written as 

(-a^+1)^^ (17) 

Here = 1 in |t), cr^ = — 1 in |^) and = 1 in |t), 
Tz = —1 in 1^). Therefore, the symmetry action on the 
boundary can be put exactly into the form of Eq. 12 with 

r(a,,a,+i) = (-<7:+i)^^^, (18) 

when the symmetry action a involves time reversal and 

r{ai,ai+i) = l, (19) 

when a does not involve time reversal. We can reorganize 
the variables and write the phase factor as a function of 

ujs{a~^ai^i,a~l^a~^,a) = /"^{ai.ai^i) (20) 

It can be checked that ujs{a~^ ai-^i^ a~^-^a~^ ^ a) is a non- 
trivial third cocycle of group Z2 x Zj, using the cocycle 
conditions introduced in appendix A. Therefore, we can 
show using methods in Ref. 3 that the boundary must be 
either gapless or symmetry breaking and the 2D bulk is 
in a nontrivial SPT phase. 



III. 3D SPT PHASE WITH Z2 x Z2 SYMMETRY 

Now we go one dimension higher and consider the 
Z2 X Z2 group. To differentiate the two Z2's, we write 
them as Z2 and Z2. We can start with similar construc- 
tions in the bulk where the nontrivial Z2 SPT phase in 
two dimension is attached to the domain wall of the Z2 
configuration in the cubes. A simple understanding of 
the ground state wave function exists starting from Levin 
and Gu's construction of 2D SPT phase with Z2 symme- 
try. It was shown that in the 2D SPT model the ground 
state wave function takes the simple form of 

i;2D{C) = (-1)^<^ (21) 

where Nc is the number of domain walls in configuration 
C. Such a wave function has gapless edge states when 
the system has a boundary. Now consider the symmetry 
Z2 X Z2. This has elements {1, ^1, ^2, ^3}- We will pick 
two of the three nontrivial elements say ^1, ^2- Choose 
any one of these two generators (say gi) and consider do- 
main walls in 3D between regions that break this symme- 
try in opposite ways. This defines closed 2D manifolds. 
Now, on these closed manifolds the domain walls of the 
second generator (^2) are examined and number of closed 
loops counted (essentially these loops are intersections of 
domain walls of gi and ^2)- Now, one uses this set of 
closed loops and defines a wave function as in Eq.21. 

V^3d(C) = {-l)^r' (22) 




FIG. 4. SPT phase in d = 3 with Z2 x Z2 symmetry Red 
(blue) surfaces represent domain walls of the Z2 (Z2) symme- 
try. They intersect along curves (black lines). The ground 
state wave function is a superposition of all domain wall con- 
figurations which differ by a sign depending on whether there 
are an even or odd number of intersection curves. This auto- 
matically implies a protected edge state when a domain wall 
intersects the surface of the sample (grey dashed line). 

This is the wave function of a Z2 x Z2 symmetric state 
where N^^^'^ is the number of closed loops formed by the 
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intersection of the gi and g2 domain walls in configura- 
tion C (see Figure 4). What is the physical consequence 
of this wave function? Consider breaking one of the Z2 
symmetries and forming a domain wall (of say ^1). Now, 
the domain wall is the edge state of the 2D SPT phase 
protected by the unbroken Z2 symmetry. This is also 
evident from the wave function. Due to this nontrivial 
topological features, the wave function describes a non- 
trivial SPT phase. 

Because there are three different ways to pick two gen- 
erators out of the three nontrivial elements in Z2 x Z2, 
we can construct three different SPT wave function in 
this way. Therefore, this construction allows us to access 
all possible nontrivial phases classified using group coho- 
mology theory. In the following sections we will present 
a more detailed study of this construction. Starting from 
a bulk Hamiltonian in section III A, we analyze its edge 
state in section III B and demonstrate its relation to non- 
trivial group cocycles. In section IIIC, we present the 
field theory description of these phases^. 



A. Ground State Wave-function on a 3D lattice 




domain wall as the original Z2 spin.^^ If two or more Z2 
domains walls touch at the Z2 vertex, then one factor is 
added for each domain wall. (We can pick a particular 
way to separate the domain walls at each vertex.) Now 
it is easy to see that the Z2 spins form 2D nontrivial Z2 
SPT states on the Z2 domain wall. Away from the do- 
main wall, they are polarized in the x direction. Now 
we include tunneling between the Z2 configurations r^^, 
which in the low energy sector of the previous Hamilto- 
nian term not only flips the Z2 spins but also changes 
the corresponding Z2 interaction pattern in its neighbor- 
hood. All these terms are local. The ground state is then 
unique and gapped and is the equal weight superposition 
of all Z2 configurations together with the Z2 SPT state 
on its domain wall and polarized Z2 spin away from its 
domain walls. The ground state is symmetric under the 
Z2 X Z2 symmetry. 



B. Connecting the edge state to Group 
Cohomology 

Imagine that we cut the system open and expose the 
boundary. Note that when doing the cut, we need to be 
careful and not to break the symmetry of the system. (In 
particular, what we can do is to double the Z2 and Z2 
spins along the cut and cut through each pair.) Then 
on the boundary, each plaquette hosts a Z2 variable 
and each vertex hosts a Z2 variable a^. It is then easy 
to see that on the ID domain wall of Z2 variables on the 
boundary is attached the ID edge state of the Z2 SPT 
phase. The Z2 symmetry acts by flipping the plaquettes 



(24) 



FIG. 5. Three dimensional SPT model with Z2 x Z2 symme- 
try. Each cube hosts a Z2 variable (big black dot) and each 
vertex hosts a Z2 variable (small green dot) (a) a cube and 
the vertices around it (b) a vertex and the cubes around it. 
Shaded surfaces represent Z2 domain walls. 



and the Z2 symmetry acts on the a variables living on 
the r domain walls as 

N 

X|cri, ...ctat) = JJ/(cri,cri+i)|o-i,...aAr) (25) 



Consider a 3D cubic lattice where each cube in the 
bulk hosts a Z2 variable (big black dot in Fig. 5) and the 
Z2 symmetry flips |0) to |1) and |1) to |0). Each vertex 
in the 3D bulk hosts a Z2 variable (small green dot in 
Fig. 5) and the Z2 symmetry flips |0) to |1) and |1) to 
|0). To define the Hamiltonian of the system, we first fix 
a Z2 configuration in the cubes and define the interaction 
between the Z2 vertices. Start with a magnetic field in 
the X direction on all the Z2 spins. If a vertex is on 
the domain wall between Z2 configurations, modify the 
Hamiltonian term ax by an extra factor of 

(_)na^^ (23) 

where Udwp is the number of Z2 domain wall pairs along 
the loop of all nearest neighbor Z2 spins on the same Z2 



where f{ai^ai-^i) adds a phase factor of i to each a do- 
main wall and is 1 if the neighboring cr's are the same 

/(ai,ai+i)=i^^^ (26) 

as given in Ref. 5. 

Equivalent ly, just like in the previous case, we can 
think of the two dimensional boundary as a system of 
independent degrees of freedom with one Z2 variable r 
per plaquette and one Z2 variable a per vertex, as shown 
in Fig. 6. The r domain walls lives on the solid lines and 
the a domain walls lines on the dashed lines. Then the 
Z2 part of the symmetry simply acts by flipping the Z2 
variables 

Z2: \ri)^\fi). (27) 
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FIG. 6. Boundary of 3D SPT phase with x symmetry. 
Each plaquette hosts a Z2 variable r and each vertex hosts 
a Z2 variable a. Circles represent local degrees of freedom 
on the 2D boundary labeled by group elements of Zi x Z2. 
Domain walls of r live on solid lines and domain walls of a 
live on dashed lines. 



The Z2 part of the symmetry action acts by flipping the 
Z2 variables gi 



Z2' \(J^)^\Gi). 



(28) 



Moreover, the Z^ symmetry also involves a phase factor i 
whenever r and a domain walls intersect. For example, at 
the shaded intersection shown in Fig. 6, the phase factor 
is given by 



I-cj-'^cjI 1- 



(29) 

Such a symmetry action can be shown to be related to 
cocycles in a similar way as in section II D. In particu- 
lar, in the group cohomology construction as reviewed in 
appendix B, the boundary local degrees of freedom are 
labeled by group elements ai of the symmetry group and 
the action of the symmetry operator involves two parts: 

• 1. changing the local states to 

• 2. multiplying a phase factor given by nontrivial 
cocycles to each triangle of a^, a/^. 

More specifically, on the 2D boundary of 3D SPT phases, 
the symmetry acts as 

0(a)|ai, ....,aAr) = '^^P{(Xi,(Xj,(X},)\(X(Xx, ...,aaAr) 

(30) 

where /^{ai^aj^ak) is a phase factor given by the non- 
trivial 4th cocycle UJ4 

^^{ai.aj.ak) = U4^{al^aj,a~^ak,a^^a~^ ,a) (31) 

To put the above symmetry action of Z2 x Z2 into 
this form, we can relate each plaquette with one of its 
vertices (the lower right corner, for example) and think 
of the Z2 and Z2 variables in this pair as labeling group 
element = (r^, a^) in the symmetry group, as indicated 
by circles in Fig. 6. Then the symmetry action labeled by 
a = T^(j flips the variables as 



Moreover, the phase factor is related to the triangle of 
ai = {Ti.Gi), aj = {Tj.Gj) and au = {rk.crk) (although 
it depends trivially on Tk and <j^). Including the depen- 
dence of the phase factor on the symmetry action applied, 
we have 

/"(ofi, Qfj, a/e) = 1 if cr G a is nontrivial 



if cr G a is trivial 



\c^i = iji.CFi)) \aai = {TTi.GCFi)) 



(32) 



Now we can reorganize the variables and check that 
a;4(a~"^aj, a) = /^{ai^aj^ak) is indeed 

a nontrivial four cocycle using the condition given by ap- 
pendix A. 

Therefore, as discussed in appendix B, the state con- 
structed in this way corresponds to a nontrivial SPT 
phase with Z2 x Z2 symmetry and the physical features 
include gapless states on Z2 domain walls on the bound- 
ary. 



C. Effective field theory description 

We access this 3+ ID topological phase by directly con- 
structing the 2 -hi dimensional boundary. While the sur- 
face is a conventional 2+ ID bosonic system except in the 
way symmetries are implemented, which prohibit a fully 
symmetric and gapped surface. We expand the physical 
symmetry Z2 x Z2 to {Z2 x Z2) x U{1). Eventually we 
will break the U{1) to Z2 and identify it with one of the 
Z2 generators. 

We begin with a bosonic field on the surface: bi ^ e^^^ 
that is charged only under the U{1) symmetry. A pos- 
sible SPT surface state is to spontaneously break this 
U{1) symmetry by condensing this boson. Restoring the 
symmetry requires condensing vortices ip of this bosonic 
field. However, for this to be the surface of a topologi- 
cal phase, the vortex transformation law must forbid a 
vortex condensate that preserves the remaining Z2 x Z2 
symmetry. This can be achieved if the transformation 
law is projective, which implies at least a two fold de- 
generacy of the vortex fields hence = (?/^+, V^-). In- 
tuitively, one may consider the elements of the group 
Z2 X Z2 = {1, ^y, ^z} as representing 180 degree 
rotations about the x, y, z axes. The projective rep- 
resentation is then just the spin 1/2 doublet, with the 
nontrivial generators represented by the Pauli matrices 
Qa = icr^ 1 thus ip icr^ip. Physical operators, which 
can be directly measured, do not realize symmetry in 
a projective fashion, however, the vortex is a non-local 
object and can hence transform projectively. Thus the 
effective Lagrangian for the vortices is: 



>Csurface = I (^M " + K{d^ay - d^a^f (33) 



where, by the usual duality, the vortex fields are cou- 
pled to a vector potential a whose flux is the number 
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density of bi bosons. Now, a vortex condensate will nec- 
essarily break symmetry since the gauge invariant com- 
bination = ip^a^ip will acquire an expectation value. 
However, this transforms like a vector under rotations 
and will necessarily break the Z2 x Z2 symmetry. Thus 
the theory passes the basic test for the surface state of 
a topological phase. This continues to be true when we 
break down the U(l) symmetry to Z2 and identify it with 
one of the existing Z2 symmetries. That is, the boson bi 
now transforms the ga- This gives us three choices, which 
can be labeled by a = X, F, Z which identifies the gen- 
erator that leaves bi invariant. Thus the phase labeled X 
has bi —bi under the generators ^y, gz but is left in- 
variant under gx - Even under this reduced symmetry the 
surface states cannot acquire a trivial gap since the con- 
densate of bi continues to break some of the symmetries. 
Now, there are three distinct nontrivial phases implied by 
this construction labeled a = X, F, Z which, combined 
with the trivial phase confirms the Z2 x Z2 classification 
of topological phases with this symmetry^ . With this re- 
duced symmetry, additional terms can be added to the 
surface Lagrangian including ACx = — (V^V^?/^ cos (/)i) 
We would like to now verify that domain walls at the 
surface carry protected modes along their length. 

Now, consider breaking the symmetry at the surface 
down to a single Z2. When the remaining Z2 generator 
is different from the one used to label the phase, we will 
see a protected mode is present along the surface domain 
walls. For example, in the setup here where the phase is 
labeled by X, consider breaking the symmetry down to 
just the Z2 generated by gz- This is realized by condens- 
ing just tp-^ or tp- vortices. Consider a domain wall in the 
x-y plane along x = 0, where for x > {x < 0) we have 
ijj-^ (V^-) condensed. The domain wall is the region of 
overlap of these condensates where we have ij\_ip- ^ e'^^^ 
which represents a vortex tunneling operator across the 
domain wall. The phase <p2 is therefore conjugate the 
boson phase (j)i as in a Luttinger liquid. Note however 
the transformation law under the remain Z2 symmetry 
is: 01^2 ^ 01,2 + TT, which is the transformation law of 
the protected edge of the 2D Z2 SPT phase^'^. A similar 
conclusion can be drawn for surface domain walls, when 
the remaining Z2 symmetry is gy- On the other hand, 
in this X phase, preserving the gx symmetry does not 
lead to protected modes since 0i is invariant under this 
generator and can be locked at a particular value without 
breaking the symmetry. 

Now, let us write down a 3+ ID topological field theory 
that describes this phase. We will restrict to an abelian 
theory, which is not ideal given that the symmetry acts 
like spin rotation along orthogonal directions. Neverthe- 
less the abelian theory gives us valuable insights. Con- 
sider the two species of bosons e*'^^ ^ introduced earlier. 
In the phase labeled X above, the number density ni con- 
jugate to 01 is invariant under all transformations, while 
the number density 77,2 changes sign under two of them. 
We can always pick these to be for group elements gx^ 
gy as above. Then we write down the following 'BF+FF' 



theory: 

27tC = e{Bidai + B2da2) + —edaida2 (34) 

27r 

where e is short for e^^^^ and the indices are sup- 
pressed. Here, 27rni = e^'^^^diBijk is the number density 
of boson 61, and Fuj = diaij—djau represents the vortex 
lines of boson bi. Thus the vector potentials ai are invari- 
ant under the transformations, but a2 changes sign under 
9x^ gv- Therefore the coefficient of the second term in 
the equation above must be quantized 6 = 0,7r by the 
usual arguments. The latter case represents the topolog- 
ical phase. The unusual properties of domain walls when 
the symmetry is broken down from Z2 x Z2 ^ Z2, at the 
surface is readily deduced from this theory. A domain 
wall which breaks gx, gv occurs when 6 = tt ^ — tt on 
the surface leads to the surface Luttinger liquid action: 

Sll ^ ^ J 9x(pidth dx dt (35) 

which describes a domain wall located along z = 0, = 
0, where we have replaced ai^2i = di<Pi,2- These are just 
the phase fields we discussed above, and in particular, 
under the remaining gz transformation they are both 
shifted by 01^2 ^ 0i,2 + tt, which is the transformation 
property of the nontrivial edge of a Z2 2D SPT phase. 

IV. GAUGING THE Z2 SYMMETRY 



Before Gauging After Gauging 




FIG. 7. Result of gauging the Z2 symmetry. The gauge flux 
now allows the domain walls to terminate on them. They 
are point particles carrying a Kramers doublet in the d = 2 
system, and flux lines that carry a one dimensional edge state 
protected by the remaining Z2 global symmetry in the d = 3 
system. 

In Ref. 5, Levin and Gu discussed the topological phase 
in d = 2 protected by just Z2 symmetry. In that case. 
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since there is no additional symmetry, we cannot inter- 
pret the domain wahs with SPT phases in d = 1. How- 
ever, there is a different mechanism whereby an SPT 
phase is generated. The domain wah ends, which are 
now particles, carry fractional statistics (semionic statis- 
tics in this case). In fact this is readily seen from the 
symmetry transformation law for the edge state of this 
phase (j)i^2 4^1,2 + TT. Now, to generate a domain wall, 
both fields must be rotated, which is achieved by the do- 
main wah operator D{x) = e^(^i(^)+^2(:c))/2^ -g readily 
seen that D{x)D{x') = zSign(x — x')D{x')D{x)^ which 
is the hallmark of semionic statistics in ID. One can 
now gauge the Z2 symmetry, to obtain a topologically 
ordered state that is distinct from the regular Z2 gauge 
theory^^'^^'^^^^. In fact the particle excitations in this 
theory are semions, and anti-semions, and this is termed 
the double-semion theory. Indeed, the act of gauging 
the Z2 symmetry simply implies the possibility of do- 
main walls ending within the 2D sample. The location of 
these ends are just the gauge charges and fiuxes. Given 
our experience with domain walls ending at the edge of 
the sample, where they have been shown to be simians, 
these end points are then expected to be semions. 

We now apply this intuition to the cased discussed 
here. First consider the d = 2 SPT phase with Z2 x 
symmetry. If the Z2 part of the symmetry is gauged, 
domain walls could end within the sample. We now ex- 
pect these ends, which are Z2 fiuxes, to carry a spin 1/2 
degrees of freedom. The gauged version this model is 
similar to the one discussed in Ref. 42. From the co- 
homology classification we know that 2D phases with 
Z2 X Z2 symmetry form a Z2 x Z2 group. Therefore, 
with this construction together with the nontrivial SPT 
phase which appears with just the Z2 symmetry, we are 
able to account for all SPT phases possible in this system. 

Now consider the d = 3 SPT phase with Z2 x Z2 sym- 
metry. Consider gauging gi the first Z2 symmetry, but 
leaving the remaining Z2 engaged so it continues to act as 
a global symmetry. Now consider and inserting a gauge 
fiux along a curve. It is readily seen that this fiux curve 
will have an SPT edge state along it, specifically, the edge 
state of the d = 2 SPT phase protected by Z2 symmetry. 
It can therefore be distinguished from a trivial phase, 
where gauging one of the Z2 degrees of freedom does not 
lead to protected excitations along Z2 fiux lines. To see 
this, note that gauging the symmetry just means that 
the domain walls of gi ends along the fiux line. However, 
this domain wall is now in a Z2 SPT phase. This follows 
from the wave function in Eqn. 22, since the loops on 
the surface bounded by this curve, formed by domains of 
^2, are weighed to give a 2D Levin-Gu wave function^. 

An interesting question is the result of gauging entirely 
the Z2 X Z2 symmetry. If we begin in the SPT phase this 
should lead to a Z2 x Z2 topological order that is dis- 
tinct from the conventional one, which described by the 
deconfined phase of a Z2 x Z2 gauge theory. Character- 
izing these subtle differences in topological order is left 
to future work. 



V. ID SPT PHASE WITH Z2 x Z2 SYMMETRY: 
FROM DECORATED DOMAIN WALLS TO 
MODEL HAMILTONIANS 



FIG. 8. An SPT phase in d = 1 protected by Z2 x Z2 symme- 
try, that act on the two sets of spins a, r. This phase emerges 
from the ordered phase of the a spins by condensing domain 
walls attached to spin flip excitations of the r spins. 

We discuss how the well known Haldane phase in ID 
can be understood within the framework of decorated 
domain walls. An advantage of this picture is that it 
leads directly to a model Hamiltonian that realizes this 
phase, and provides a simple rationale for the string order 
parameter of this state. Consider breaking down the full 
SO (3) spin rotation symmetry to just Z2 x Z2 symmetry, 
which is sufficient to define this topological phase. 

The two Z2S can be modeled by a pair of Ising mod- 
els cr, r, and the ordered phases are given by (a^) 7^ 
or (r^) 7^ or {a^r^) 7^ 0. Consider beginning in the 
partially ordered phase (a^) 7^ but (r^) = 0. Now, 
the domain walls of a and the Z2 quanta of r are both 
gapped. If we condense the former, we enter the com- 
pletely disordered phase. Condensing the latter leads to 
the completely ordered phase. But one can consider the 
following scenario. What happens when we condense the 
bound state of a domain wall and the r Z2 quanta (see 
figure 8)? 

It is readily seen that this describes the SPT phase. 
First, since we condense domain walls of a, the Z2 sym- 
metry is restored. Note, although r quanta are con- 
densed, they are condensed along with the domain walls, 
so this is not a local operator, which implies that Z2 is 
also unbroken. So we have the full Z2 x Z2 symmetry. 
The easiest way to see this is the SPT phase is to con- 
sider the 'order' parameter for this phase, which is just 
the condensate p = a^r^, where is the disorder param- 
eter that creates/destroys a domain wall in the a fields. 
So we expect long range order in: 



n '^"(^) 

i<k<j 



However this is just the string order parameter for the 
Z2 X Z2 SPT phase"^^. 
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The following space-time picture motivates why it has 
edge states. Consider the path integral representation in 
imaginary time, with a spatial boundary. Then space- 
time is a cylinder, the periodic direction being time. The 
condensate of Z2 'charged' domain walls leads to world 
lines that sometimes intersect the boundary. When they 
do, since they carry Z2 quantum number, they flip the 
spin at the edge - so there must be gapless edge degrees 
of freedom that fluctuate in time. 

This picture also motivates the following Hamiltonian. 
First consider a pair of decoupled Ising Models, a and r 
as shown in the figure. To enforce the binding of charge 
to domain walls we add the following Hamiltonian: 

H, = -XJ2 + TiLl/2TiVl/2<) (36) 

i 

The Hamiltonian J^puster44,45 j^g^. ^ g^^. q£ commut- 
ing projectors, and has a unique ground state on a sys- 
tem with periodic boundary conditions. However, gapless 
edge states appear when the system is terminated at a 
boundary. 

VI. 3D SPT PHASE WITH TIME REVERSAL 
SYMMETRY FROM DECORATED DOMAIN 
WALLS PICTURE 

The decorating domain wall picture allows us to con- 
struct nontrivial SPT phases with time reversal sym- 
metry also. In this section we are going to discuss 3D 
phases with symmetry only, x U{1) symmetry 
and Z2 X 5*0 (3) symmetry. In particular, the phase we 
construct with symmetry only is beyond the coho- 
mology classification. 

Consider a 3D lattice with a spin 1/2 sitting in each 
cube. Time reversal symmetry maps between spin up 
It) to spin down \D states together with the complex 
conjugation operation in this basis. Now we can decorate 
the domain wall in this spin configuration with some 2D 
states. First we need to specify the orientation of the 
domain walls as pointing from \D to Then we can put 
chiral states on the domain wall whose chirality matches 
the orientation of the domain walls. In particular, we 
put the Kitaev's Eg state^^ on the domain wall which is 
a bosonic 2D state with no fractional excitations in the 
bulk and a c_ = 8 edge state. After this decoration, we 
sum over all possible spin configurations, together with 
the ^8 decoration. This wave function is invariant under 
time reversal symmetry which can be seen as follows. 
Time reversal maps between ||) and |t) changing the 
orientation not the potision of the domain walls. At the 
same time, time reversal maps the state to —^8 state 
with inverse chirality. Therefore, the chirality of the Eg is 
always consistent with the orientation of the domain wall 
and the total wave function is time reversal invariant. 

Such a construction gives rise to a 3D SPT state with 
only time reversal symmetry. Following similar argu- 



ments as in previous sections, we see that on the bound- 
ary of the system if we break time reversal symmetry in 
opposite ways on neighboring regions, then the domain 
wall between these regions carry a chiral edge state with 
c_ = 8. This is half of what one can get in a pure 2D sys- 
tem without fractionalization. This is consistent with the 
field theory construction with BE + EE term in Ref.9. 
On the other hand, if time reversal symmetry is not bro- 
ken, then the surface is gapless or has topological order. 
In fact, the surface topological state is that with three 
fermions^^ which cannot be realized in pure 2D system 
with time reversal symmetry. 

Similar construction applies to 3D SPT phases with 
X U{1) and Zj x S0{3) symmetry. In 2D with U{1) 
and SO (3) symmetry, there is an integer class of SPT 
phases with even integer quantized charge or spin Hall 
conductance.^' ^^'^^'^^ Similar to the construction above, 
we can put the first nontrivial phase in this class to the 
domain wall of time reversal. The Hall conductance chi- 
rality should match the domain wall orientation. By sum- 
ming over all spin configurations, we obtain a state with 
both time reversal and U{1) or 5*0 (3) symmetry. One 
signature of the state would be even integer quantized 
Hall conductance on the time reversal domain wall on 
the 2D surface of the system. 

VII. RELATION TO KUNNETH FORMULA 
FOR GROUP COHOMOLOGY OF GROUP 

G = GixG2 

If a group G is the direct product of two subgroups 
G = Gi X G2, then the Kiinneth formula for the group 
cohomology of G can be written as^^ 

H^+^(G, U{1)) = H^(Gi, H^+^-^(G2, U{1))) (37) 

k=0 

which says that cohomology groups of G in (i + 1 di- 
mension can be obtained from cohomology groups of Gi 
and G2 in lower dimensions. Due to the relation be- 
tween cohomology groups and SPT phases, this implies 
that SPT phases with G symmetry in d dimension can be 
constructed from SPT phases with Gi and G2 symmetry 
in lower dimensions. 

When /c = 0, the term in the formula is 

n'>{Gi,n''+HG2,u{m = n''+HG2,u{i)) (38) 

which means that some SPT phases with symmetry G in 
d dimension are just SPT phases with symmetry G2 in d 
dimension. 

When /c = d + 1, the term in the formula is 

n''+\Gun°{G2,u{m = n''+HGuU{i)) (39) 

which means that some SPT phases with symmetry G in 
d dimension are just SPT phases with symmetry Gi in d 
dimension. 
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Our domain wall construction correspond to the term 
with k = 1 

H\Gi,n\G2^U(l))) (40) 

1-L^{G2-,U{1)) labels SPT phases with G2 symmetry in 
d - 1 dimensions. If H^(G2,/7(1)) = M (M = Z^, Z 
for example), then the term becomes l-L^{Gi^ M) which 
labels one dimensional representations of Gi using M 
coefficient. Therefore, this term says that some d dimen- 
sional SPT phase with G symmetry can be obtained from 
d—1 dimensional SPT phases with G2 symmetry and one 
dimensional representation of Gi with M coefficient. 

Our domain wall construction provides a physical in- 
terpretation of the above statement. With discrete Gi 
group, as in the cases we are interested in, consider a 
d dimensional configuration with group elements in Gi. 
The domain walls are then also labeled by group elements 
91 ^ Gi given by the difference of the group elements on 
the two sides of the wall. Then on the d—1 dimensional 
domain wall labeled by gi , we can choose to put d—1 
dimensional SPT phases from the class M . Our choice 
should be consistent with the fusion rules of the domain 
walls. That is, the phases and G M which we 
choose to put onto domain walls and g\ should com- 
bine into the phase m^^ which we put onto the domain 
wall labeled by ^i^^. In other words, this correspond to a 
one dimensional representation of group Gi in coefficient 
M. 

Therefore, the term with A: = 1 in the formula 37 cor- 
respond exactly to our domain wall construction. Note 
that this formula also gives the condition of when such 
domain wall construction can be consistent. In particu- 
lar, when putting the d—1 dimensional SPTs onto the 
domain walls, we need to put them according to the cor- 
responding one dimensional representation of Gi in M. 
Suppose we have M = Z3, then putting the ffist nontriv- 
ial one onto a Z2 domain wall is not allowed because two 
Z2 domain walls can merge to trivial and so should the 
corresponding SPTs. 

Terms with higher /c's would correspond to construct- 
ing d dimensional SPT phases with G symmetry by 
putting d — k dimensional SPT phases with G2 symmetry 
onto d — k dimensional defects in Gi configurations. 



VIII. CONCLUSION AND DISCUSSION 

In the previous sections we have presented the con- 
struction of a two dimensional SPT phase with Z2 x 
symmetry by attaching Haldane chains, a one dimen- 
sional SPT phase with Zj symmetry, to the Z2 domain 
walls in the 2D bulk and also a three dimensional SPT 
phase with Z2 x Z2 symmetry by attaching the 2D SPT 
state with Z2 symmetry to the domain wall of the Z2 vari- 
ables in the 3D bulk. Such a construction leads directly 
to the special topological feature of Z2 domain walls on 
the boundary of the system: the Z2 domain walls on 



the boundary carry gapless edge states of the other sym- 
metry {Z2 and Z2). We established the nontrivial SPT 
order in the system by relating the symmetry action on 
the boundary of the system to nontrivial cocycles and 
also demonstrated how the SPT order can be properly 
described using field theories. 

This domain wall construction also applies to 3D SPT 
phases with Zj, Zj x U{1) and Zj x 50(3) symmetry, 
including the time reversal symmetric topological 'super- 
conductor' phase, which features chiral edge modes 
along surface domain walls that break time reversal sym- 
metry. This state lies beyond the group cohomology clas- 
sification. The domain wall construction of attaching 
d—1 dimensional SPT states of G2 symmetry io d — 1 
dimensional defects in Gi configurations was shown to 
be related to one term in the Kiinneth formula for the 
group cohomology of groups of the form Gi x G2. Other 
terms in the formula may be related to attaching lower 
dimensional SPT states with G2 symmetry to lower di- 
mensional defects in Gi configurations, exploring which 
is left to future work. A ID SPT phase with Z2 x Z2 
symmetry, when approached in this manner, readily sug- 
gests a parent Hamiltonian. It would be interesting to 
find physically well motivated Hamiltonians that realize 
the higher dimensional topological phases as well. The 
physical viewpoint on SPT phases described in this work 
may help guide such a search. 
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Appendix A: Projective representation and group 
cohomology 

Matrices u{g) form a projective representation of sym- 
metry group G if 

u{gi)u{g2) uj{gi,g2)u{gig2), ^1,^2 ^ G. (Al) 

Here uo{gi^g2) ^ ^(1) and oj{gi^g2) 7^ 1, which is called 
the factor system of the projective representation. The 
factor system satisfies 

^^^^'H^2, ^3)^(^1, ^2^3) = uj{gi,g2)uj{gig2,g^), (A2) 

for all ^1,^2,^3 ^ G^ where 5(^1) = 1 if ^1 is unitary and 
5(^1) = -1 if ^1 is anti-unitary. If ^2) = 1, V^i,^2, 
this reduces to the usual linear representation of G. 
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A different choice of pre-factor for the representation 
matrices u\g) = /3{g)u{g) wih lead to a different factor 
system u\gi,g2): 



P{9i)P'^''H92) 
P{9i92) 



^{91^92) 



(A3) 



We regard u\g) and u{g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems oo'{gi^g2) and oj{gi^g2) as be- 
longing to the same class uj. 

Suppose that we have one projective representation 
ui{g) with factor system uJi{gi^g2) of class uji and an- 
other U2{g) with factor system 002(91^92) of class ob- 
viously ui{g)<^U2{g) is a projective presentation with fac- 
tor system cji (^1, ^2)^2(5'i, ^2)- The corresponding class 
(J can be written as a sum uji-\-uj2. Under such an addi- 
tion rule, the equivalence classes of factor systems form 
an Abelian group, which is called the second cohomology 
group of G and is denoted as 1-C^[G^ The identity 

element 1 G V?[G^ is the class that corresponds to 

the linear representation of the group. 

The above discussion on the factor system of a projec- 
tive representation can be generalized which gives rise to 
a cohomology theory of groups. 

For a group G, let M be a G-module, which is an 
abelian group (with multiplication operation) on which G 
acts compatibly with the multiplication operation {ie the 
abelian group structure) on M: 

g • {ah) = {g • a){g -6), g e G, a,b e M. (A4) 

For the cases studied in this paper, M is simply the U{1) 
group and a a /7(1) phase. The multiplication operation 
ab is the usual multiplication of the U{1) phases. The 
group action is trivial g-a = a{geG,ae U{1)) if g is 
unitary and ^ • a = a* if ^ is anti-unitary. 

Let ujn{gi gn) be a function of n group elements 
whose value is in the G-module M. In other words, 
ujn ' G"" ^ M. Let C^[G,M] = {un} be the space of 
all such functions. Note that C^[G^M] is an Abelian 
group under the function multiplication uj'^{gi^ ...^ g^) = 
uJn{9ii"'i9n)^'n{9ii"'i9n)' We define a map dn from 
C^[G, /7(1)] to C^+i[G, /7(1)]: 

((in^n)(^l, ...,^n+l) = 

9l ' ^n(^2, ^n+l)^i"^^"^' (^1, "',9n) X 
n 

'[\'^'k^^\9i,-,9i-i,9i9i+i,9i+2,-gn+\) (A5) 



1=1 



Let 



S"[G,M] = {w„K = d„. 



(A6) 



and 



S^[G,M] and 2:^[G,M] are also Abelian groups which 
satisfy S^[G,M] c Z^[G,M] where B\G,M] = {1}. 
Z^[G, M] is the group of n-cocycles and S"^[G, M] is the 
group of n-coboundaries. The nth cohomology group of 
G is defined as 

[G, M] = [G, M] /S^ [G, M] ( A8) 

When n = 1, we find that uji{g) satisfies 

ooii9iK^''\92)^oJi{9i92) (A9) 

Therefore, the 1st cocycles of a group are the one dimen- 
sional representations of the group. 

Moreover, we can check that the consistency and equiv- 
alence conditions (Eq. A2 and A3) of factor systems 
of projective representations are exactly the cocycle and 
coboundary conditions of 2nd cohomology group. There- 
fore, 2nd cocycles of a group are the factor systems of the 
projective representations of the group. 

When n = 3, from 



{dsUJ3){9i,92,9s,94) 
_ ^s{92, 93^ 94)^s{9i^ 929s, 94)^3(91^ 92, 9s) 
^3 {9192 ,93,94)^3 {91 1 92 , 9394 ) 

we see that 
Z^[G, [/(!)] = {W3 1 

A91) I 



(AlO) 



(All) 



^3 {92,93, 94)^3 {91 ,9293, 9a)^3 {91 ,92,93) ^ 

^3 {9192 ,93,9^)^3 {91 ,92,9394) 



and 



B'[G,U{l)]={co3\co3{9i,92,93) 



<^2^'\92,93)<^2i9l,9293) 

^2(5152, 53)^2(51, 52) 

(A12) 



Z^[G,M] = {LUn\dniUn = l,u;„ S C"[G,M]} (A7) 



which give us the third cohomology group H^IG, U{1)] = 
Z3[G,/7(1)]/B3[G, [/(!)]. 



Appendix B: Group cohomology and symmetry 
action on the edge state of SPT phases 

In establishing the nontrivial SPT order for the mod- 
els discussed in section II and III, we used the connection 
between group cocycles and the symmetry action on the 
edge degrees of freedom. Here we review this connection 
established in 1 where the specific form of the edge sym- 
metry action is given in terms of cocycles for the exactly 
solvable models. Because of such an nontrivial form of 
symmetry action, the edge state of nontrivial SPT phases 
cannot be trivially gapped without breaking symmetry. 
Note that in generic models, symmetry action on the 
edge degree of freedom can take more general form and 
it is not completely clear (in 3D and higher) how to ex- 
tract the cocycle information from the symmetry action. 
However for the models discussed in section II and HI, 
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the edge symmetry action takes exactly the form given 
in 1 and therefore allows us to identify the cocycle and 
establish the nontrivial SPT order in a straight forward 
way. 

ID SPT phases are classified by projective represen- 
tations of the symmetry group G with factor system 
^2(^1,^2) ^ [/(!)]. If the edge degree of freedom 

has basis states labeled by group elements |^), g G G, 
then one possible form of the projective representation is 



U{g)\gi 



UJ2 



{9i ^9 \9)\99i) 



(Bl) 



That is, the symmetry action changes group element ba- 
sis from gi to ggi and also adds a nontrivial phase fac- 
tor given by uj2{gi^ g~^ g) • It can be checked that U{g) 
does form projective representation with factor system 
^2(^1,^2). 

2D SPT phases are classified by the third cocycles of 
the symmetry group G, 6^3(^1,^2,^3) eU^[G,U{l)]. The 
edge of the system can be thought of as a local ID system 
with a special symmetry action. If the degrees of free- 
dom on the ID edge have basis states labeled by group 
elements gi G G, then one possible form of the effec- 
tive symmetry action on the edge is 



U{g)\gi,g2--^gN) 



\99i^992,'",99n) 



9i+i9~ 



That is, the symmetry action changes group element ba- 
sis from gi to ggi and also adds a nontrivial phase factor 
between every pair of nearest neighbor states given by 

^3(^^^^^+l'^i+l^"^^)• 

3D SPT phases are classified by the fourth cocycles of 
the symmetry group G, uj4{9ii 92, 93^ 94.) ^ n^[G,U{l)]. 
The edge of the system can be thought of as a local 2D 
system with a special symmetry action. To describe the 
symmetry action, we need to triangulate the surface such 
that the degrees of freedom live on the vertices of the 
triangulation. If the degrees of freedom have basis states 
labeled by group elements gi G G, then one possible 
form of the effective symmetry action on the edge is 

U{g)\gi,g2---,9N) YlA^^i97^9j,9j^9k,9k^9~^,9)x 
\99i,992,'",99n) 

That is, the symmetry action changes group ele- 
ment basis from gi to ggi and also adds a nontrivial 
phase factor on each triangle on the surface given by 
^A{9i^9j,9j^9^,9k^9~^,9)' 

Symmetry action on the edge taking the particular 
forms given above is a sufficient but not necessary condi- 
tion for the existence of SPT order. In general, the edge 
degrees of freedom do not have to have group element ba- 
sis and the symmetry action can take more general form. 
However, for the discussion in this paper, these particu- 
lar forms are enough to establish the SPT order in our 
models. 



Appendix C: Symmetry action on low energy states 
of Eq.5 from exact diagonalization 

In section II B, we studied one particular realization 
of the edge dynamics on the boundary of the two- 
dimensional SPT phase with Z2 x Z2 symmetry. The 
Hamiltonian governing the edge dynamics is given by 



He 



and the Z2 symmetry acts on the edge as 

i 

while the symmetry acts as 



(Cl) 



(C2) 



-i(jl)K (C3) 



In this section, we show how to extract the low energy 
effective action of the symmetry on the edge state from 
exact diagonalization. 

The Hamiltonian is an XY model (written in xz plane 
here) on a spin 1/2 chain and the low energy effective 
theory is known to be described by the compactified free 
boson theory with Lagrangian 



2llCedge = St0l9x02 - V 



(9x02) 



(C4) 



The low energy states are labeled by quantum numbers 
Uk ^ Z and n/e G Z, /c = 0, 1, 2..., where no and no label 
the total angular momentum and the winding number of 
the boson field respectively and and n/e, /c > 0, label 
the left/right moving boson modes. If we normalize the 
ground state energy to be and the first excited state 
energy to be 1/4, then the energy of each low energy 
state is given by 



E = 



k>0 



k{nk +n/c) 



(C5) 



and the lattice momentum p {—ir/a < p < n/a) of each 
state is given by 



2k, 



k>0 



(C6) 



where L is the total system size, no is given by the con- 
served U{1) quantum number 



(C7) 



The boson field (j) = ^^'^^'^ can be thought of as de- 
scribing the direction the spin 1/2 is pointing to in the 
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xz plane. Therefore, it corresponds to the spin state 
cos^lt) + sin 011). From this, it is easy to see that the 
Z2 symmetry maps (j) to —cj). However, it is not easy to 
see the action of the time reversal symmetry on the low 
energy state because of its complicated form. In order 
to obtain this information, we perform exact diagonal- 
ization of the Hamiltonian, identify the n/c, fik quantum 
numbers of the eigenstates from quantities like p, Sy 
and find out how time reversal symmetry acts on these 
states. 

Energy levels in the spectrum can be degenerate and 
p and Sy allows us to partially split the degeneracy. 
However, states with {no = 0, no, {n/c, n/c}} and {no = 
0, — no, {^/c, ^fe}} have the same p, and Sy. In order 
to tell them apart, we need to use our knowledge of the 
action of the Z2 symmetry. Using the decomposition of 
the 6 field 



{x,t) = 



TTQVt 

L 



L 



E 



he 



-ip{x-vt) _^ ^] ^ip{x-vt) _^ ^^^ip{x+vt) _^ ^] ^-ip{x+vt) 



where tto measures total angular momentum (with 
quantum number no), tto measures winding number 
(with quantum number no) and hk are the anni- 
hilation operators for left/right moving boson modes 
(with occupation number labeled by Uk and n/c). Be- 
cause the Z2 symmetry maps (j) to —0, it maps tto, 
TTo, ^k{^\)i ^/c(^I) all to minus themselves. There- 
fore, under Z2 symmetry action, |no, no, n/c, n/c) goes to 
(_)Efc>o |— 77,Q^ — no, {n/e, n/e}). Therefore using the 
action of the Z2 symmetry, we can further distinguish 
states with {no = 0, ±no, {^/c, ^/c}} and fix the relative 
phase between them. If we want to fix the global phase 
factor of these two states, we can calculate the action of 
complex conjugation on them, which is expected to act 
as |no,no,{n/c,n/c}) to |-no, no, {n/e, n/e}). 

Now we are ready to look at each degenerate sec- 
tor labeled by p, Sy and see how symmetry acts on 
them. We discuss the first few sectors as an illustration 
of method. The spectrum is obtained by diagonalizing a 
system with 16 spin 1/2's. 

First the ground state with = 0, p = 0, = is 
nondegenerate. If the phase factor is fixed such that the 
state is unchanged under complex conjugation, then it is 
invariant under both Z2 and time reversal symmetry. 

The first excited states with E = 1/4 and p = 
are two-fold degenerate, one with Sy = 1 and one 
with Sy = —1. The two states are |1,0,{0,0}) and 
| — 1,0, {0,0}). If we fix gauge such that complex conju- 
gation acts as in this subspace and Z2 symmetry 



acts as 



1 

1 



in this subspace, then time reversal acts 



-1.0000 -0.0065i 0.0000 - O.OOOOi \ 
1^ 0.0000 - O.OOOOz -1.0000 - 0.0095Z ) ' 
The third energy level at p = is also two fold 
degenerate with states |2,0,{0,0}) and |-2, 0, {0, 0}). 



If we fix gauge such that complex conjugation acts 
as in this subspace and Z2 symmetry acts as 

in this subspace, then time reversal acts as 

0.9619 - 0.0004i 0.0381 - 0.0099A 
0.0381 - 0.0099Z 0.9619 - 0.0004z ) ' 
The third energy level also has two fold degeneracy 
at p = 1 with states |0, 1,{0,0}) and |0, -1, {0, 0}). 
If we fix gauge such that complex conjugation acts 

as in this subspace and Z2 symmetry acts as 

in this subspace, then time reversal acts as 

-0.0380 + 0.6342Z -0.9621 + 0.0251A 
-0.9621 + 0.025H -0.0380 + 0.6342z ) ' 
The third energy level also has two fold degeneracy at 
p = '^^j2 with states |0,0, {ni = 1, others are 0}) and 
|0,0, {ni = 1, others are 0}). If we fix gauge such that 

complex conjugation acts as in this subspace and 



Z2 symmetry acts as 



time reversal acts as 



in this subspace, then 



-1 
-1. 
-1' 
-1 

Therefore, up to finite size inaccuracy, the action of 
time reversal on the low energy state is consistent with 
|no,no,{n/e,n/e}) to (-)^'«>o |no, -no, {n/e, n/e}). 



Appendix D: Chern-Simons theory description of 2D 
SPT phases 

In this section we briefiy review the Chern-Simons ap- 
proach to 2+1-D SPT phases introduced in Ref. 6. In 
particular the U(l) x U(l) Chern-Simons effective the- 
ory and effective edge theory of bosonic Z2 x Zj-SPT 
phases will be discussed. 

It is believed that any Abelian gapped phases (all quasi 
particle excitations obey Abelian statistics) in 2+1-D 
can be described by U{1)^ Chern-Simons theory^^^^. 
In general the bulk effective theory for a 2+1-D gapped 
Abelian phase is 



C = 



47r 



(Dl) 



where K is a real symmetric N x N matrix of integer 
arguments and • • • represents non-universal higher-order 
terms (such as Maxwell terms). The ground state de- 
generacy of the gapped phase is | detK|^ on a Riemann 
surface of genus g. For bosonic SPT phases which only 
support bosonic excitations in their spectra (no aynons), 
we must require detK = ±1. On the other hand, the 
gauge invariance of the bulk Chern-Simons theory (Dl) 
on an open manifold suggests the existence of edge exci- 
tations on the manifold's boundary^^: 



(D2) 
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where V is a positive-definite real symmetric matrix 
which depends on microscopic details of the system. {^/} 
are chiral boson fields which is quantized to have 27r pe- 
riodicity: 0/ + 27r. The chirality of edge modes 
is determined by the signature (n+,n_) of matrix K: 
each positive (negative) eigenvalue of K corresponds to a 
right (left) -mover on the edge. Hence for the edge states 
to be non-chiral with the same number of right- and left- 
movers, n+ = n_ and dimK = N =even are required. 
Besides, for a bosonic system the diagonal elements of 
matrix K are all even integers while for a fermionic sys- 
tem, at least one diagonal element of K are odd integers. 

A key ingredient for the Chern- Simons approach to 
SPT phases is the implementation of symmetry in ef- 
fective theory (Dl) and (D2). For a non-chiral gapped 
phase, in the absence of symmetry its edge states in 
general are always gapped due to backscattering. Such 
backscattering terms on the edge has the following 
generic form of 

ACedge = ^1 COs(^ + Ci). (D3) 

1 / 

where 1 = {h^hr ' ' ^^n) are integer vectors and C\ 
are all real parameters. These backscattering terms cor- 
respond to Higgs terms which condense bosonic quasi 
particles^^ in the bulk^. In the presence of symmetry, 
however, certain backscattering terms (D3) can be for- 
bidden due to nontrivial symmetry transformations of 
chiral bosons {(/>/}. In other words, symmetry can pro- 
tect the edge states from backscattering, yielding gapless 
edge excitations. Since the universal properties of SPT 
phases are fully encoded in their boundary excitations, 
the fingerprints of a SPT phase is contained in the ma- 
trix K and how the edge chiral bosons transform under 
symmetry. In general for a symmetry group G, the chiral 
bosons transforms under symmetry operation ^ G G as 

g^ig-' = v^g ^ j Wf >j + Sc^j, (D4) 
(W^)^KW^ = r]9K, 0<S(l)i < 27r. 

where 77^ = ±1 for unitary (anti-unitary) symmetries and 
G GL{N, Z) is di N X N unimodular matrix. These 
symmetry transformations {W^, 5(j)^} on the edge chiral 
bosons form a faithful representation of symmetry group 
G for the 2+1-D SPT phase. And the fuh set of data 
[K, {W^, characterizes a SPT phase with symme- 

try group G. 

With the same K matrix, different symmetry transfor- 
mations sometimes correspond to the same SPT phase. 
How to identify different SPT phases in terms of their 
data [K, {W^, It turns out with preserved sym- 

metry group G, different gapped phases always form an 
Abelian group. The identity element of this group is 
the trivial phase, whose edge excitations can be gapped 
out without breaking any symmetry. For example, the 
state whose edge chiral bosons don't transform under 
symmetry belongs to the trivial phase. The other ele- 
ments of the group are different SPT phases. The natu- 
ral Abelian product rule of two phases [K, {W^, Scj)^}] 



and [K, {W^, 6(j)^}] is to take their matrix direct sum: 

[K,{w^(5(/)n]e[K,{w^^^n] 

= [K © K, {W^ © © 6^9}], (D5) 

It's easy to figure out the inverse of a phase is 
[K,{W^(^(^^}]-i = [-K,{W9^S(j)9}]. 

Therefore two data [K, {W^ and [K,{W^^0^}] 

will represent the same phase if and only if 

[K, {W^ S(l)9}] 0[-K, {W^ S(j)9}] = trivial. (D6) 

In the following we illustrate the above general princi- 
ples in an example: bosonic SPT phases with Z2 x 
symmetry. We denote the generator of Z2 and time re- 
versal symmetry as X and T respectively. For a bosonic 
SPT phase, the simplest K matrix is the following 2x2 
one 

K^(;;). ,0.) 

It turns out this 2 x 2 K is enough to describe all different 
bosonic SPT phases with Z2 x Zj symmetry, which form 
an Abelian Z2 x Z2 group. These SPT phases transform 
in the following way under Z2 x symmetry^: 

T ( T-1 = ^ -^1 ^ 

\h J \h + riTTT J 

where nx = 0, 1 and tit = 0, 1 correspond to the 4 differ- 
ent Z2 X Zj^-symmetric gapped phases. Let's label these 
4 phases as [nx,^T]- Among them [0,0] is the trivial 
phase, whose edge states can be gapped out by backscat- 
tering term A cos((/)2 + G). All the other 3 phases are non- 
trivial SPT phases, whose edge modes cannot be gapped 
out without breaking Z2 x symmetry. However, it's 
straightforward to show that 

[nx^nr] ©[nx, ^t] = trivial. 

and hence these 4 phases form a Z2 x Z2 group. 

In the following we prove that the phase [0, 1] with 
symmetry transformations (E7) has a different but equiv- 
alent representation (E5)-(E6). Consider a [0,1] phase 
with edge chiral bosons {^1,2} transforming as (E7), is 
put together with a phase ^ whose edge chiral bosons 
{^1,2} transform as (E5)-(E6). It's easy to check that the 
following backscattering terms which preserves Z2 x 
symmetry can be added to the edge: 

Ai sin(^i +^1), 
Ai sin(^i - 0i), 

A2 cos (j)2 COS (j)2- 

The first two terms localize the value of chiral boson fields 
01 and (j)i up to tt, while the last term fixes this ambigu- 
ity. When these three backscattering terms all become 
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relevant, the edge states will be fully gapped out with- 
out breaking any symmetry. This means [0, 1] i> = 
[0, 1]0[0, 1] =trivial and hence <l = [0,1]. And we've 
shown that (E5)-(E6) is just a different but equivalent 
representation of the SPT edge with symmetry transfor- 
mations (E7). 



From these low-energy operator correspondence and sym- 
metry operation (2)- (3), it's clear the chiral bosons trans- 
form as 



X 



(E5) 



Appendix E: Field Theory for the Edge Lattice 
Model of a 2D SPT Phase 



under Z2 operation X and 



In this section we derive the low-energy long- 
wavelength effective theory for the edge states of Hamil- 
tonian (4). Since this model describes the 1+1-D edge 
states of a gapped 2-hl-D state, periodic boundary con- 
dition is always imposed. Therefore the usual duality 
transformation^^ needs to be modified, by introducing 
an extra Z2 gauge field^ /i^. To be specific, we perform 
the following duality transformation to the spin- 1/2 vari- 
ables and a* in 1+1-D Hamiltonian (4): 



,,z \x \x 

P2i-l,2i^2i-1^2i5 



M2i-l,2i5 



^x — P2i,2i+1^2i^2i+l5 — P2i,2i+1' 



A? 



-A? 



(El) 



^2i-l5 'z^z — ^^2i- 

And the gauge- invariant constraints are 

^2n-lM2n-2,2n-lM2n-l,2n — ^5 
^2nf^2n-l,2nf^2n,2n-\-l = fl- 
it's straightforward to see that ni(~l)*'^f — ^- Since 
there are L (length of the 1+1-d system) independent 
constraints on dual variables jli-i^i and A^, these dual 
variables faithfully represent the 2^-dimensional Hilbert 
space of original spin- 1/2 variables and <j\ And the 
Hamiltonian (4) can be rewritten by dual variables as 



(E2) 



The low-energy effective Lagrangian of such a XY model 
can be easily obtained^'^^ as 

27r£erfge = dt4>idA2 - v[K{^f + l(a,(/.2)2](E3) 

where Luttinger parameter K = 1 m model (4) and v = 
4a is the velocity of edge excitations (a denotes the lattice 
constant). In the low-energy the chiral boson fields {^1,2} 
are related to local spin operators by 



TT 2a 2a ' 



,i(/)2(a?) 



it:' 



(E4) 



' M2n-l,2n'^2n-1^2n ^ ^2n ,2n+1^2n^2n+l 



- T«(l - a 



cr. 



,n+l 



X X ) 



f TT 
- TT 



(E6) 



under time reversal operation T. As shown in Appendix 
D, such a symmetry transformation for the 1+1-D edge 
state of a 2 + 1-D Z2 x Zj-SPT phase has another equiv- 
alent but different representation: 



X 
T 



(t>2 
62 + TT 



(E7) 



In the new representation (E7) it's easy to figure out the 
following operator D creates a Z2 domain wall: 



D{x) ^ e 



i02(a^)/2 



(E8) 



since we have the following commutator 



[^i(x),e^^^2^^^] = i/7r-sgn(^-x)-e^^^2(^\ W u e R. 



due to the Kac-Moody algebra [(/)i(x), 9^02(^)1 = 
27ri(5(x — y) from Lagrangian (E3). It's straightforward 
to show that 



Hence when time reversal acting twice the domain wall 

operator obtains a minus sign: D[x) —D{x). There- 
fore the low-energy effective theory of the edge states also 
predicts that a Z2 domain wall on the edge of Z2 x Zj- 
SPT phase should carry the projective representation of 
time reversal symmetry, just as implied in the lattice 
model (4) with symmetry (2)- (3). 
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